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INTRODUCTION 


The approach of the project is to base the design of multi-function, reflective topographies 
on the theory that topographically dependent phenomena react with surfaces and 
interfaces at certain scales. The first phase of the project emphasizes the development of 
methods for understanding the sizes of topographic features which influence reflectivity. 
Subsequent phases, if necessary, will address the scales of interaction for adhesion and 
manufacturing processes. 

A simulation of the interaction of electromagnetic radiation, or light, with a reflective 
surface is performed using specialized software. Reflectivity of the surface as a function 
of scale is evaluated and the results from the simulation are compared with reflectivity 
measurements made on multi-function, reflective surfaces. 


METHODS 

Simulation 

In this work a numerical simulation of light interaction with a surface is compared with 
reflectivity measurements made at NASA with a Perkin-EImer Lambda5 scatterometer. 
Light is emitted from a source, reflected by a surface and then intersects a detector. 

Topography 

In the simulation, we represented the topography with a 200 x 200 grid of points acquired 
by a scanning tunneling microscope (STM), where the points are located in a grid in x and 
y with a height z. The STM work was performed at NASA Langley on a Digital 
Equipment Nanoscope II. Six topographies were scanned at three scan sizes, 20jj.m x 
20pm, 2pm x 2pm, and 200nm x 200nm (see Table 1). 

The large, topographic data sets, used to represent the surface, are analyzed by the 
patchwork method where the surface, represented by the data points, is tiled with 
triangular patches (Brown et al. 1992). The topography is evaluated over a range of 
scales, or patch sizes, by tiling over the surface with decreasing patch sizes. In the 
simulation each patch represents a reflective facet, atomically smooth and a perfect 
reflector. 

The triangular patches are placed on the surface in two directions: parallel and 
perpendicular to the STM scan direction. Reflectivity calculations are made for each 
direction and the results of the simulation are the average of the two calculations. 

Lambda5 Scatterometer 

We created a computer model of the Perldn-Elmer LambdaS scatterometer from optical 
path representations of the scatterometer and reflectivity assembly which were provided 
by Perkin-EImer. The incident angle of the light is user-defined, and the detector was 
modeled as a rectangle at the position and orientation defined by the optical schematic. 
The size of detector was defined from an engineering drawing of the detector which was 


provided by Hammamatsu (part number R298 HA). The output of the detector was 
assigned a value of unity for any intersection of a reflected ray with the detector. 

Table 1 STM Scan and Reflectivity Measurements 


.FACE 

STM FILENAME 

MATERIAL 

POINT SPACING 

REFLECT % 




X x Y x Z (nm) 

20 DEGREES 





531 nm 

A 

NASA22 

T7M1R-P1 

100 

69.1 

A 

NASA23 

T7M1R-P1 

10 

69.1 

A 

NASA24 

T7M1R-P1 

1 

69.1 

B 

NASA25 

P4H1R 

100 

46.7 

B 

NASA26 

P4H1R 

10 

46.7 

B 

NASA27 

P4H1R 

1 

46.7 

C 

NASA16 

D25M2R-P30 

100 

44.7 

C 

NASA17 

D25M2R-P30 

10 

44.7 

C 

NASA18 

D25M2R-P30 

1 

44.7 

D 

NASA19 

T7H3R 

100 

44.1 

D 

NASA20 

T7H3R 

10 

44.1 

D 

NASA21 

T7H3R 

1 

44.1 

E 

NASA 13 

D25H1R-P30 

100 

15.2 

E 

NASA14 

D25H1R-P30 

10 

15.2 

E 

NASA15 

D25H1R-P30 

1 

15.2 

F 

NASA10 

P4M1R 

100 

14.1 

F 

NASA11 

P4M1R 

10 

14.1 

F 

NASA12 

P4M1R 

1 

14.1 


The light emitted by the source was modeled as set of parallel rays that originate from the 
source, travel to the surface and are reflected by the center of each patch. One ray is 
generated for each patch. The direction of the reflected ray is calculated from the incident 
ray direction and the normal of the patch. From the center of the patch, the ray is 
reflected off of an optical wedge, a concave, spherical mirror and then it is determined if 
the ray intersects the detector. 

The output of the Lambda5 is percent reflection relative to a known reference sample, in 
this case the reference sample was a stainless steel mirror provided by Perkin-Elmer. The 
reflectivity measurements are expressed in terms of a percentage of the measured 
reflectance from the reference sample. 

Reflectivity Simulation 

We ran the computerized simulation on the 18 STM data sets with an incident angle of 20 
degrees. The simulation generated the incident rays, reflected them off the patches and 
counted the rays which rays intersected the detector. The output of the simulation R, or 
absolute percent reflectivity, is defined as the number of incident rays that intersect the 
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detector is divided by the total number of rays reflected by the surface. R is calculated for 
each patch size and is plotted versus log(patch area). R r , or relative reflectivity, is 
calculated from R as 


Rra-b (RA-RB)/RA [eq. 1] 


where RA>RB. 

The simulation was modified to account for the effect of patch orientation on the intensity 
of the incident rays. The intensity of the incident ray, initially equal to one, is multiplied by 
the cosine of the angle between the patch normal and the incident ray. The ray is then 
reflected by the patch and the simulation records the number of rays and their intensities 
that intersect the detctor. The output of the weighted reflectivity simulation, R w is 

R w = (Z intensity collected rays) / number of reflected rays. [eq. 2] 


Reflectivity Results 

We combined the simulation's reflectivity results from the three scans sizes into one larger 
set of results using Matlab, a matrix-based software program. The large set of results, 
combined for each of the six surfaces, covered a range of patch areas of 7 orders of 
magnitude, from 0.5 nm 2 to 3pm 2 . Table 2 shows the range of scale for each of the three 
scan sizes and the data point spacing of the scans. 

Table 2. Large and Small Patch Sizes for Three STM Scan Sizes 


Scan Size 
20um x 20um 
2 um x 2 um 
200 nm x200nm 


Large Patch Size Small Patch Size 
3 pm 2 5000nm 2 

5000nm 2 50 nm 2 

50 nm 2 0.5 nm 2 


Point Spacing 
lOOnm 
lOnm 
lnm 


When the scan sizes were combined, the magnitudes of R did not correspond at the joining 
patch size. The values of R were shifted by the difference between the scans so as to 
match at the joining patch size. The R values of the largest scan size (20um x 20um) were 
used as the zero shift scan when the scans were shifted up, and the R values of the two 
scans (2pm x 2pm and 200nm x 200nm) were shifted up to the zero shift scan. When the 
scans were shifted down, the smallest scan size (200nm x 200nm) was the zero shift scan, 
and the two larger scans were shifted down to it. The shifted results were used to 
generate plots of absolute reflectivity vs. log(patch area). 

Relative reflectivity results were calculated from the shifted scans for shifted up results, 
and plots of relative reflectivity vs. log(patch area) were generated. Calculation of relative 
reflectivity (Rg), from eq. 1, was designed to factor out the dependence on the reference 
sample. Surface A, which has the largest reflectivity measurement, was used as the 
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reference in the plots. A negative value of R r indicates that the reflectivity of the surface 
is greater than surface A. 

Scale of Interaction 

The results of the reflectivity simulation are compared to experimental results, obtained 
from NASA Langley's Lambda5, to calculate a scale of interaction of the light with the 
surface. The scale of interaction was defined as the square root of the patch area, from 
the reflectivity simulation, where the corresponding magnitude of R is equal to the 
reflectivity value measured by NASA on the Lambda5. Figure 1 shows how a scale of 
interaction is found from the simulation and experimental results. 


REFLECTIVITY (ABSOLUTE) vs SCALE 



LOG (Scale of Measurement 
or Triangle Patch Size) 


Figure 1- Sche matic of Calculation of Scale of Interaction from Idealized, Absolute 
Reflectivity vs. Scale The intersection of the measured reflectivity is found, and the scale 
of interaction is calculated from the corresponding patch size. The scale of interaction is 
shown to occur in region 2. 


RESULTS 
Relative Area 

The STM data sets were analyzed by the patchwork method and a representative scale- 
area plot is shown in figure 2 (Brown et al. 1993). The scale-area analyses were 
conducted on each of the three scan sizes separately. All of the scale-area plots are found 
in Appendix A. 
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Absolute Reflectivity 

Absolute reflectivity as a function of scale, or patch area, was calculated by the reflectivity 
simulation. Shifted up an shifted down results for surface A are shown in Figures 3 and 4. 
Absolute reflectivity plots for the six surfaces are found in Appendix B. 

Cosine Weighting - incident angle 

Absolute, cosine-weighted reflectivity as a function of scale, for surface B, is plotted with 
non-weighted results in Figure 5. Cosine-weighted results for the six surfaces are found in 
Appendix C. 

Relative Reflectivity 

Relative reflectivity as a function of scale, for surface B relative to A, is shown in figure 5. 
Relative reflectivity plots for surfaces B through F, relative to A, were generated using 
shifted up results. The plots are found in Appendix C. 

Scale of Interaction 

A scale of interaction was calculated from the absolute reflectivity plots. Table 3 lists the 
scales of interaction found using the absolute reflectivity results. The reflectivity 
measurements did not intersect the relative reflectivity results, and no scales of interaction 
were found. 


\.data\l angleysnasa 1 1 Quadrant <1,1) of <1,1) 

200 relative areas calculated Crossover lou: 205 

Aver agerPatch Size a Surface Coverage Crossover high: 2.8ie+03 

Relative Surface Area vs. Patch Area Data-dependant 



Figure 2 - Scale-Area Plot of Surface A. The 2mm x 2mm scan size is shown for 
surface A. The relative area begins to increase (crossover) at a patch area of 28 10 nml 
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REFLECTIVITV (ABSOLUTE) vs SCALE - SHIFTED UP 



Figure 3 - Shifted Up Absolute Reflectivity Results for Surface A The percent 
reflectivity is 100% at large patch areas and decreases to about 53% at the fine patch 
sizes. The largest scan size was used as the no shift scan. 


REFLECT IWITY (ABSOLUTE) vs SCALE - SHIFTED DOUH 



Figure 4 - Shifted Down Absolute Reflectivity Results for Surface A. The percent 
reflectivity is about 48% at large patch areas and decreases to about 0% at the fine patch 
sizes. The smallest scan size was used as the no shift scan. 
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REFLECTIVITY (ABSOLUTE) v« SCALE - SHIFTED UP 



Figure 5 - Incident Angle Cosine- Weighted and Non-Weighted Absolute Reflectivity 
Results for Surface B. The difference between the cosine weighted and unweighted 
results is shown. Cosine weighting shifts the set of results down by a factor of 0.94 
(a multiplier of cosine 20°), or shifting the results down by 6%. The results are shifted by 
an equal percentage at all scales. 


REFLECTIVITY (RELATIVE) vs SCALE - SHIFTED UP 



Figure 6 - Relative Reflectivity Results for Surface B Relative to Surface A. The 
maximum relative reflectivity value is about .15 which is less than the measured result. R r 
is negative at patch areas less than about 500 nm 2 - 
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surface with a scale of interaction from the shifted up results. The scales from the shifted 


down differ by about 1 order of magnitude. 


Surface 

Scale of Interaction (nm) 

Scale of Interaction (nm) 


shifted up 

shifted down 

A 

39 

x 

B 

X 

X 

C 

X 

102 

D 

X 

164 

E 

X 

13 

F 

X 

18 


DISCUSSION 
Relative Area 

The scale-area plots show that the reflective topographies are complex at fine scales, and 
the relative areas, a measure of complexity, increase with decreasing patch size. The 
difference between the relative areas of the three scan sizes is clear: the maximum, 
relative areas of the 200nm x 200nm scan size are approximately 2 orders of magnitude 
larger than that of the 2 pm x 2pm scan size and approximately three orders of magnitude 
larger the smallest scan size. 

Absolute Reflectivity 

The plots show three distinct regions over a range of scale. In region 1, occurring at 
patch sizes down to 10 5 nm 2 , the percentage reflectivity (R) remains constant over a range 
of patch sizes. Decreasing patch size does not change the amount of reflected light that 
reaches the detector, and R is largest in this region. 

In region 2, occurring at patch sizes from 10 5 nm 2 to 100 nm 2 , R decreases with 
decreasing patch size. Decreasing patch size decreases the amount of reflected light that 
reaches the detector, and R in this region is less than region 1 and greater than region 3. 

In region 3, occurring from a patch size of 100 nm 2 to 0.5 nm 2 , R remains constant. 
Decreasing patch size no longer decreases the amount of light that reaches the detector, 
and R is lowest in this region. 

The material properties of the surfaces, i.e., conductivity, absorbtivity and transmissivity, 
were not considered in the simulation, and would shift the plots down from the shifted up 
condition. The downward shift would increase the calculated scale of interaction, and, 
depending on the size of the shift, would cause the measured reflectivity results to 
intersect the simulation's results. 
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Incident angle cosine weighting shifts the reflectivity simulation results down for all of the 
surfaces, and the percentage shift is constant for all scales. It is speculated that the 
percentage shift is constant at all scales because the collector is small, and only rays close 
to the direction of specular reflection are collected. Reflected rays with a weighting 
factor close to zero would reduce the results by more than 6%, but these rays, with a small 
weighting factor, will not be reflected close to the specular direction, and will not intersect 
the detector. The weighting factors of the collected rays are close to cosine (20°), and 
shift the results down by 6 %. 

Relative Reflectivity 

The measured results do not intersect the relative reflectivity results at any scale, which 
may be because the scale of interaction theory is wrong or because we are misinterpreting 
the results of the reflectivity measurements. It was expected that all relative values would 
be positive, at least at one scale, because the calculations were made relative to the surface 
with the largest reflectivity measurement, surface A, and that this scale would correlate 
with the reflectivity measurements. The plots show that surface A is less reflective than 
most of the other surfaces and that the surface with the largest R r changes with scale. 

The reflectivity measurements made by NASA are expressed as a percentage reflectivity of 
a stainless steel reference sample. It is not yet clear how we should interpret this 
representation of reflectivity compared to the computer simulation. Relative reflectivity 
plots were generated to factor out the dependence on the reference sample, but we have 
not been successful. Including conductivity and absorbtivity of the reference sample in the 
computer simulation may provide a truer representation of NASA's reflectivity method. 
Also, more information about how the LambdaS processes the output signal from the 
photo multiplier tube may give a better understanding of the equipment's output. 

Scale of Interaction 

Joining the three scan sizes effects the scale of interaction calculation; shifting up 
decreases the scale and shifting down increases it. Possible causes of mismatch at joining 
patch areas are the patch placement algorithm, differences in the STM scan parameters for 
the three scan sizes or variability in material properties over the different scan areas. 

The current algorithm places a small number of patches at the large scales (large with 
respect to data point spacing), and may provide a poor representation of the topography. 
Since the topography may be more precisely represented with decreasing patch size, the 
joining patch size may be thought of as a boundary separating regions of high and low 
precision. 

Changes in material properties of a surface will also change how the tip interacts with the 
topography. Efforts were made to minimize the effect of local changes in material 
properties, but the scans may have been effected to some degree. 
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Future Work 

The work will be continued under a NASA training grant, grant number NGT-5 1 107. In 
future work other reflectivity methods, such as total integrated scattering (TIS), will be 
investigated as a means for better understanding the amount of energy reflected by the 
surfaces. Reflectivity samples will made from homogeneous materials to reduce the 
complexity of the reflectivity simulation by eliminating multiple layer materials and 
distributed reflective particles. Random patch placement algorithms will be investigated 
that may better represent the interaction of light with the surfaces. 

CONCLUSIONS 

1. The scale-area plots show that the reflective topographies are complex at fine scales, 
and that the relative areas, a measure of complexity, increase with decreasing patch 
size. 

2. Simulation of a reflecting surface as a collection of triangular mirrors and decreasing 
the size on each repetition results in a steady decrease in the amount of light arriving 
at a detector, indicating increasing scatter or diffuse reflection at finer scales. 

3. Ranking of the surfaces based on reflectivity calculated from the current algorithm 
does not correspond at any scale, to the ranking from reflectivity measurement, as 
they are currently interpreted. 

4. The scales of interaction calculated from the current algorithm do not share a 
common region of reflectivity with the measured values. 
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Appendix A 
Scale-Area plots 




\data\langlei 
200 relative 
Average/Patcl 
Rel 


1.00 uZ 

r z =0. 863423 
g =-0.031859 
6. 186e+01 to 




\data\langlei 
200 relative 
Average/Patcl 
Rel 


r z =0. 737951 
y =-0. 159053x+0. 63953 
1.297e+03 to L297e+0* 




Relative Area 


\dat a\l angl eg\nasa 1 3 

200 relative areas calculated 

Average/Patch Size & Surface Cover; 

Relative Surface Area vs. 


1.001 



1-000 t : , 

10,000 

r *=0.628521 

y =-0. 000398x+0. 00 1 966 
7. 136e+03 to 7. 136e+04 







\dat a\lang 
200 relati 
Average/Pc 

1.9 
1.8 
1.7 
1.6 
1.5 
1.4 
1.3 
1.2 
1.1 


1.0 

r *=0.99855 
y =-0.1009 
5.076e-01 
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\data\lang 
200 relati 
Average/Pc 


r 2 =0. 84693 
g =-0.069-1 
6.098e+01 



\data\langlei 
200 relative 
Average/Patcl 
Rel 


r 2 =0.995190 
y =-0.224827 
2.499e+0G to 










Relative Area 








\dat a\l angl eg\nasa22 

200 relative areas calculated 

Aver age/Patch Size & Surface Coverage 

Relative Surface Area vs. Patch Area 


Quadrant (1,1) of (1,1) 
Crossover lou: 7.02e+03 
Crossover high: 1.19e+05 
Dat a-dependant 



5.316e+03 to 5.316e+04 



Relative Area 


\data\l angle 
200 relative 
Average/Pati 
Re 


1.0 


r 2 =0. 891796 

y =-0.051 28 lx+0. 169147 
6.096e+01 to 6.096e+02 


1,000 

Patch 


10,000 

Area Cnm 2 ) 




\data\lang 
200 relati 
Average/P< 
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r z =0.9921 3 
y =-0.331! 
1.890e+01 



\data\langley\nasa25 

200 relative areas calculated 

Average/Patch Size & Surface Coverage 

Relative Surface Area vs. Patch Area 


Quadrant (1,1) of (1,1) 
Crossover lou: 1.32e+04 
Crossover high: 2.46e+05 
Data-dependant 
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\data\langl< 
200 relativi 
Aver age/Pat 


r z =0. 960182 

y =-0. H7334x+0.707816 
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Appendix B 

Absolute Reflectivity vs. Patch Area 
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Appendix C 

Cosine Weighted and Non-Weighted Absolute Reflectivity vs Patch Area 
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Appendix D 

Relative Reflectivity vs. Patch Area 
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